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Abstract. We study the evolution of growth and decay laws for the magnetic field coherence length ^, energy Em and 
magnetic helicity H in freely decaying 3D MHD turbulence. We show that with certain assumptions, self-similarity of the 
magnetic power spectrum alone implies that ^ ~ t^^^. This in turn implies that magnetic helicity decays as ~ t^^" , 
where s — (^diff/C-tf)^, in terms of ^diff, the diffusion length scale, and ^h, a length scale defined from the helicity power 
spectrum. The relative magnetic helicity remains constant, implying that the magnetic energy decays as Em ~ ^-1/2-2^ jjjg 
parameter s is inversely proportional to the magnetic Reynolds number Rem, which is constant in the self-similar regime. 
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1. Introduction 

Magnetic fields are ubiquitous in the Universe, being ob- 
served in objects from planets to galaxy clusters (Zeldovich 
et al. 1983, Ruzmaikin et al. 1988, Kronberg 1994). In galax- 
ies and galaxy clusters, the typical strength is of order a few 
/iGauss, which is thought to be produced by dynamo action 
on a seed field. In galaxies the dynamo timescale is roughly a 
rotation period, 10^ yr, and a simple calculation (Ruzmaikin 
et al. 1988) based on the age of a typical galaxy shows that 
the seed field must have been about lO^'^^ Gauss, or perhaps 
less in the currently favored models with a cosmological term 
(Davis et al. 1999). 

There is no shortage of ideas for generating this seed 
field. The more conventional astrophysical explanations are 
based on a Biermann battery operating at the era of reioniza- 
tion (see, e.g., Gnedin et al. 2000, and references therein). 
There are more speculative ideas based on various gener- 
ation mechanisms in the early Universe (Grasso & Rubin- 
stein 2001), which have the common feature of producing 
stochastic, homogeneous and isotropic magnetic and velocity 
fields, characterized by their power spectra and initial length 
scales. Some of these mechanisms produce stochastic fields 
with non-zero magnetic helicity (Joyce 1997, Cornwall 1997, 
Vachaspati 2001) the first of these being maximally helical. 
Another common feature of these generation mechanisms in 
the early Universe is that they last for a short time only, typi- 
cally much less than the time it takes for the Universe to dou- 



ble in size at the time of generation, after which the magnetic 
fields decay. As we are discussing mechanisms operating at 
the era of the electroweak phase transition (10^^^ s) or be- 
fore, field is generated essentially instantaneously compared 
with any current astrophysical or cosmological timescale. 

The subsequent decay of these primordial fields, and 
also of those in star-forming regions, motivates the study 
of freely decaying magnetohydrodynamic (MHD) turbulence 
(Mac Low et al. 1998, Biskamp & Muller 1999, Muller & 
Biskamp 2000, Christensson et al. 2001). The decay will not 
in general just be through linear dissipation, as the fields in 
the early are likely to have high magnetic Reynolds number 
because of the high conductivity of a fully ionized relativis- 
tic plasma (we recall that Rem = £,v/ri, where ^ and v are 
the typical length scale and velocity of the system, and rj the 
conductivity). This has been taken to mean in the cosmolog- 
ical context that the magnetic field is frozen into the plasma, 
and the scale length of the field increases only with the ex- 
pansion of the Universe. This is in general untrue, because 
the plasma can move, and turbulence can transfer energy to 
different length scales (Brandenburg et al. 1996). 

Thus a study of decaying MHD turbulence is required in 
order to calculate quantities such as the size of the seed for the 
galactic dynamo or the amplitude of the perturbations in the 
temperature of the Cosmic Microwave Background (CMB) 
radiation arising from primordial magnetic field generation 
(Durrer et al. 1998, 2003, Caprini & Durrer 2001, Caprini et 
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2 MHD EQUATIONS 



al. 2003). Our numerical results (Christensson et al. 2001) un- 
covered decay laws for magnetic fields which were not those 
one would expect from purely linear dissipative processes. 
In particular we saw that helical fields decayed more slowly 
than non-helical, due to the fact that magnetic helicity is an 
invariant of ideal MHD. Helicity is known to be important in 
dynamo theory (Pouquet et al. 1976, Meneguzzi et al. 1981, 
Brandenburg 2001), and we shall also be able to confirm its 
importance in decaying turbulence. Our interest here is to try 
and understand the results, and to compare them with those 
found by Biskamp & Miiller (1999) and Miiller & Biskamp 

(2000) . In doing so we have developed a new framework for 
understanding scaling in decaying 3D MHD turbulence, in 
the case where the fields are close to being maximally heli- 
cal, as in the mechanisms proposed by Joyce & Shaposhnikov 
(1997) and Vachaspati (2001). It should be noted that none of 
the estimates by Durrer et al. (1998, 2003), Caprini & Durrer 

(2001) , and Caprini et al. (2003) take into account the decay 
laws we find, and so our results are of direct importance for 
cosmology. 

The decay of the magnetic field in the turbulent case is of- 
ten presumed to result from an inverse cascade (Pouquet et al. 
1976, Meneguzzi et al. 1981, Brandenburg 2001), in which 
power is transferred locally in fc-space from small to large 
scales. However, while it is certainly true that power is trans- 
ferred from small to large scales, as can been seen from the 
energy power spectra plotted in Christensson et al. (2001), 
it was not established that the power is transferred locally, 
despite the appearance of the term "inverse cascade" in the 
title of our paper, and so it may be strictly incorrect to call 
the process a cascade. In fact, a true inverse cascade seems 
rather unlikely, as large scale power appears almost immedi- 
ately even from initial conditions which are highly localized 
in fc-space. For our analysis in this work it will not matter 
whether or not there is a cascade, and we will not discuss 
the matter further. However, we emphasize that the decay of 
the energy is not simple linear dissipation: there is definitely 
interesting non-linear physics, as shown by the appreciable 
Reynolds numbers (100-600) and the processing of the initial 
power spectra. 

Various scaling arguments have been put forward to ob- 
tain the growth law for the length scale of the magnetic field 
and the decay law for the energy. For ideal MHD (infinite 
conductivity), Olesen (1997) and later Son (1999), Field & 
Carroll (2000), and Shiromizu (1998) argued ^{t) - t2/("+5) 
where t is conformal time, and n is the index of the ini- 
tial magnetic power spectrum. Supporting evidence for this 
scaling law was also found in two-dimensional MHD simu- 
lations (Galtier et al. 1997). Shiromizu's results were based 
on a renormalization group argument, which has been revis- 
ited by Berera & Hochberg (2001), who do not find evidence 
for an inverse cascade in the steady state. The effect of having 
significant helicity was supposed to modify this scaling law 
to ~ t^/^ Em t-^^'^ (Biskamp 1993, Son 1999, Field 
& Carroll 2000). Early numerical experiments with a shell 
model of the full MHD equations (Brandenburg et al. 1996) 
suggested ^ t^-^^, and gave supporting evidence to the in- 
verse cascade. Full MHD simulations by Biskamp & Miiller 



(1999) and Muller & Biskamp (2000) showed, in the heli- 
cal case, an energy decay law Em ~ t^^^^, supported by a 
phenomenological model, which we will discuss at the end of 
this work. Decaying non-helical turbulence was studied semi- 
analytically in a shell model by Basu (2000) who found an 
energy decay law of Em ~ The decay of helical fields 

was also studied semi-analytically by Sigl (2002), giving a 
growth law for the length scale between t^/^ and t^^^. Un- 
fortunately, direct comparison with our results is otherwise 
difficult because the correlation functions are expressed in 
real space. The importance of magnetic helicity in slowing 
down the decay has been recognized earlier in studies of a 
low-order model of three-dimensional hydromagnetic flows 
(Stribling & Matthaeus 1991), where it was also found that a 
finite initial cross helicity (not studied in the present paper) 
can slow down the decay. However, power law exponents of 
the decay have not been determined for this low-order model. 

In an earlier paper (Christensson et al. 2001) we per- 
formed 3D simulations both with and without magnetic he- 
licity, starting from homogeneous and isotropic random ini- 
tial conditions, with power spectra suggested by cosmolog- 
ical applications. We found that the coherence scale of the 
field grows approximately as t^^^, with significant transfer of 
power to small scales in the helical case, which we ascribed to 
an inverse cascade. The magnetic power spectrum was self- 
similar with an approximately fc^'^-^ behavior at high fc. We 
found decay laws for the magnetic and kinetic energies of 
i^"^ and t^^-^ in the helical case, and t~^-^ for both in the 
non-helical case. These are close to, but not identical to those 
found by Biskamp & Muller (1999), and we suggested that 
their relatively large initial length scale, 25% of the simula- 
tion box size, might account for the difference. It should be 
emphasized that we are interested in the decay only while the 
scale length of the flows is less than the simulation volume, 
as we want results relevant to the early Universe where there 
are no boundaries. 

In this paper we present a new theoretical understanding 
of our numerical results for the power law behavior of the 
length scale and the energies in the helical case. We show that 
the key to understanding the power laws is the self-similarity 
of the magnetic field, coupled with the near-invariance of the 
helicity, and that the crucial parameter controlling the rate of 
decay of the magnetic energy and the helicity is the magnetic 
Reynolds number 

Our theoretical model has analogies with decaying fluid 
turbulence in 2 dimensions, where there is also an ideal in- 
variant, the kinetic energy, which plays a similar role to the 
helicity in 3D. Decaying turbulence in 2D (Ting et al. 1986, 
Chasnov 1997) exhibits self-similarity, and power-law decays 
in the kinetic energy and enstrophy (mean squared vorticity) 
are observed in numerical simulations at high Reynolds num- 
ber (Chasnov 1997). We reserve detailed discussion for Sec- 
tion|5l 

2. MHD equations 

The matter and radiation in the early Universe is modeled as 
an isothermal compressible gas with a magnetic field, which 
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3 3D MHD SIMULATIONS OP DECAYING TURB ULENCE 



is governed by the momentum equation, the continuity equa- 
tion, and the induction equation, written here in the form 
(Waleffe 1993) 

— = -u- Vu - dVhipH 

dt p 



dlnp 
dt 
dA 
'dt 



— u • V Inp — V • u, 
u X B + ryV^A, 



(1) 
(2) 
(3) 



where B = V x A is the magnetic field in terms of the mag- 
netic vector potential A, u is the velocity, J is the current 
density, p is the density, p is the dynamical viscosity, and 
T] is the magnetic diffusivity. In an expanding Universe the 
equations are identical when expressed in terms of confor- 
mally rescaled fields B, u and dissipation parameters i^, and 
rj. We work in the gauge A" = —rj'V-A (Subramanian & 
Barrowl998, Christensson et al. 2001) under the assumption 
that T] is uniform. 

In Fourier space it is useful to represent the vector poten- 
tial in terms of its projection onto an orthogonal basis formed 
by e+, e_ and k, or 



Ak = A+e+ + A e_ + Alk. 



(4) 



The two basis vectors and e_ can be chosen to be the 
unit vectors for circular polarization, right-handed and left- 
handed respectively. That is e± = ei ± ie2 where ei and 62 
are unit vectors orthogonal to each other and to k. They are 
given by ei = kxz/|kxz| and §2 = kx (kxz)/|kx (kxz)| 
respectively, z is a reference direction. 
Note that since 



(5) 



where s = ±1, this corresponds to an expansion of the mag- 
netic vector potential into helical modes. 

Using these basis vectors it is easily seen that the shell- 
averaged magnetic energy spectrum is 

EMik) = 2ne{\B^\^) (6) 



where the amplitude of the magnetic field is given by 

iBkp^dA+p+iA^nikp 



(7) 



and the expression for the shell-averaged magnetic helicity 
spectrum H{k) is 

=47rfc2(A* -Bk) (8) 

where 

A^Bk = (|A+p-|A^|2)|k|. (9) 
The function H{k) is bounded in magnitude by the inequality 

\H{k)\ <2k-^EM(k). (10) 

A field which saturates the above inequality is maximally he- 
lical. 



3. 3D MHD simulations of decaying turbulence 

We solve Eqs. ([l)-(|3} numerically (Brandenburg 2001) using 
a variable third order Runge-Kutta timestep and sixth order 
explicit centered space derivatives. All runs are performed on 
a 120^ grid, using periodic boundary conditions, as is appro- 
priate when modeling an infinite volume system. The average 
density (p) = po (where the brackets denote a volume aver- 
age) is conserved. 

We use natural units where the speed of light is c = 1, 
and fix the unit of length by setting ki — 1, where fci is 
the smallest wave number in the simulation box. Hence the 
box has size 27r. The scale factor is fixed by setting po = 1, 
and B is measured in units of ^popoc, where /io is the mag- 
netic permeability. We define the mean kinematic viscosity as 
V = pI Po. The sound speed Cs = l/\/3, as appropriate for a 
relativistic fluid. 

The equations are not quite those for a relativistic gas in 
the early universe (Brandenburg et al. 1996). However, we 
have checked that our results change little when using the true 
relativistic equations in the low velocity limit (Christensson 
etal. 2001). 

We take u and B to be homogeneous and isotropic Gaus- 
sian random fields drawn from a power-law distribution with 
a high wavenumber cut-off. The mechanisms for the produc- 
tion of magnetic fields, such as the helical production mech- 
anism of Joyce & Shaposhnikov (1997), are all stochastic, 
homogeneous and isotropic, and have an associated length 
scale k^^. Hence the power spectra, Pmik) 
and Py{k) = (u^ • Uk) can be initially modeled as 



(Bj; Bk) 



(11) 

Py{k) = Ay/^"^e-('=/'==)^ (12) 

Causality demands that um > 2 and ny > (Durrer et 
al. 1998, 2003, Caprini & Durrer 2001, Caprini et al. 2003). 
Note that in the plots it is the shell-integrated energy spectra, 
EM,v{k) — Ank^ x ^PM,v{k), which are shown. 

The amplitudes A^ can be chosen independently, pro- 
vided A*_^ = A^, which is just the condition that the vector 
potential be real. Therefore it is possible to adjust the ampli- 
tudes \A^ I and \A^ \ freely and in so doing obtaining a mag- 
netic field with arbitrary magnetic helicity. With our method 
we are able to put statistically random but maximally helical 
fields in our initial conditions. In our runs with initial helicity 
we take H — iJmax- We took tim.v to have the lowest values 
consistent with causahty, and chose kc ~ 30. 

The initial magnetic energy was taken equal to the kinetic 
energy, and had the value 5 x 10^"^ in all runs. The initial 
density was uniform and equal to 1 . The values for v and ?/ 
for the runs presented here are summarized in Tabled 

These values were chosen so as to maximize the Reynolds 
numbers while maintaining numerical stability and resolv- 
ing the dissipation scales 2'K\lvi and 2Tr^/r]t. The magnetic 
dissipation scale is resolved for times t > l/N'^r], where 
N = 120 is the lattice size, with a similar formula for the 
viscous damping scale. We start taking data at i ~ 40, when 
^diff — 27r\/8 for the smallest of our resistivity values, which 
corresponds to fc ~ 30, or approximately 3 lattice points. 
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4 HELICITY CONSERVATION AND MAGNETIC ENERGY DECAY 



RunABCDEFG 4 



Wu To OTf 07 OTf OTf OTf 07 

lO^V 1-0 0.7 0.6 0.5 0.4 0.3 0.2 

vt/^H 0.49 0.59 0.62 0.64 0.66 0.67 0.69 

lO'^ReZ/ 9.2 5.4 4.5 3.7 3.0 2.4 1.6 
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2.5 
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1 

0.5 

^0 30 40 50 60 70 80 

t 

Fig. 1. Time evolution of the ratio of length scale to the dif- 
fusion scale for run D. The notation is i^i/^diff (dashed), 
^H/Cdiff (dot-dashed), £,R/^dis (dotted) and ^xZ-^diff (contin- 
uous). The length scales are defined in Eq. J14> and the sub- 
sequent text. 

the helical length scale goes as ~ t^^"^, the decay law for 
the magnetic energy is 

Eu-^t-^/^-^'. (18) 
In any case, the energy cannot decay faster than this. Given 
^^REL — £,r/£,i, it is seen from Fig.^that Hy\el is indeed of 
order unity and does not decay markedly with time. 

To characterize the decay laws we define the exponents 

Qit) = -tEu/EM, Rit) = -tH/2H. (19) 
In Fig. |2]we have plotted R{t) versus the quantity s{t) = 
(Cdiff/Cff for several runs with different initial conditions. 
The time span is t ~ 40 to ~ 80, i.e. the last half of 
the simulations. This figure tells us several things. Firstly, it 
shows us that the value of R is approximately independent 
of time which confirms the power law decay of H. Secondly, 
Fig. 121 indicates that the quantity s is also approximately in- 
dependent of time, hence reinforcing the relation ^ t^^^- 
Thirdly, it seems that R and s are almost equal, and not just 
proportional. 

Taking the relative helicity to be constant, it follows from 
the power-law behavior of H that the energy decay law is 
indeed Em ~ t^^/^^^'*. In the limit of exact conservation of 
magnetic helicity, s ^ 0, the magnetic energy must decay as 

Em ~ t-^/'. 

Regarding the physical significance of the parameter s we 
note that if — vt, where v is the RMS velocity, (i.e. if 
the eddy turn-over time is t) then s ~ (27r)^/i?eM, where 
ReM is the magnetic Reynolds number evaluated using the 
helicity scale ^h- We have measured / = ft/^H and Re^^ 
for all runs, and find that they both change by less than about 
10% between < = 20 and the end of the runs at t = 80, 
giving the final values in Tabled One can see that there is a 
linear relation / = /o + /i / Rem between the two, and a least 
squares fit gives /o = 0.734 ± 0.002, /i = -26.6 ± 0.4. 



H = ~2V I ^,kmAi\'-\A^\'). (13) 



Table 1. Kinematic viscosity jy, magnetic diffusivity 77, 
vt/S^n and the inverse of the magnetic Reynolds number 
^eivi = ^nv/r] for our runs, is the helicity length scale 
defined in Eq. \\A\ and v is the RMS velocity, evaluated at 
the end of the runs {t ~ 80). 



Examining Fig. 1 of Christensson et al. (2001), we can see 
the magnetic energy spectrum turning down for fc ^ 30 at 
t = 46.3, showing that the system is exhibiting Ohmic dissi- 
pation as required. Under-resolved runs show ringing in real 
space and the velocity and the magnetic field explode fairly 
quickly, and so are in practice easy to discard. 

4. Helicity conservation and magnetic energy 
decay 

In this section we shall take a closer look at the observed scal- 
ing laws for helical magnetic fields and see if we can under- 
stand them on theoretical grounds. Magnetic helicity is not 
exactly conserved if 77 7^ 0, as it straightforward to show that 

(2^)3 

Hence, providing the power spectrum of the gauge field de- 
cays faster than for k ^ 00, implying that the magnetic 
energy spectrum Euik) decays faster than k^^, we can de- 
fine a helicity scale such that 

ij = -8Ai//e|. (14) 

If we assume that the evolution of is described by a power 
law ^ it is clear that the solutions to Eq. ilAl are qual- 
itatively different depending on the exponent r. If and only if 
r = 1/2 does the magnetic heUcity show a power law decay 

H^t-'^" (15) 

where 

s = (Cdiff/Cff)', (16) 
in terms of the diffusion scale ^diff = 27rY^. Addi- 
tional length scales we consider are the integral scale = 
27T J dkk^^EM{k)/ J dkEuik), the relative helicity scale 
^R, — tt\H\/Em and the magnetic Taylor microscale = 
2TTBrms/ Jims, where Srms and Jmis are the RMS magnetic 
field and current density respectively. It is plausible that all 
these scales are proportionally related and Fig. [l] shows that 
this is indeed the case. A theoretical reason for this behavior 
is given below. 

The magnetic helicity is bounded in magnitude by the in- 
equality Eq. ilO\ . An equivalent way to express this is by 

_tt\Hm[ 

-HreL = -r-^^ — < 1, (1/) 

where iJREL is the relative magnetic helicity. So if the above 
bound remains approximately saturated, i.e. i?REL ^ 1, and 
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Fig. 2. The quantity R{t), defined by Eq. il9\ versus the 
quantity s{t) — (^diff/C-ff)^ for runs listed in Table [2 The 
length scales and ^diff are defined in Eq. il4\ and the sub- 
sequent text. The time span is approximately from i ~ 40 to 
~ 80, i.e. the last half of the simulations. 
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Fig. 3. Energy decay exponent Q, defined in Eq. ( I19l l. plotted 
against the inverse magnetic Reynolds number Re^, calcu- 
lated using the helicity scale ^h, for the last half of the runs 
listed in Tabled The fit is a quadratic least squares fit using 
the final values for each run. 



must scale as t^^^, assuming the scaling form of the magnetic 
energy power spectrum 

£^M(fc) = rVf(fcO (20) 

found in Christensson et al. (2001). From the MHD equations 
one can show 



Em = - (fxu - {3 xB) -7] / d^xJ 



(21) 



We assume that the dissipation term, if not dominant, always 
contributes a constant fraction to the energy loss. If Em (k) ^ 
k^^ at high k, with z < 3, then the integral is dominated by 
the high-fc cut-off ■^J;^. We do indeed see a power spectrum 
with 2 < z < 3 (Christensson et al. 2001), which means 
that there is a well-defined helicity decay law, and that the 
energy decay is dominated by the smallest scales. Hence the 
dissipative loss e^/ — —Em is then 



CM 



3+z 



diff 



(22) 



However, on integrating the energy power spectrum and then 
differentiating with respect to time, 

eM-iC'-^ -t-^C"-'- (23) 

If we use the fact that ^diff ^ t^^"^, we find on equating the eM 
found in the two different ways, ~ i-(i+2)/2^-?-i^ 

or 



1/2 



(24) 



In deriving the power law for the length scale of the magnetic 
field, it is not necessary that the dissipative terms account for 
all the energy loss, just that it scales the same way as the 
energy loss and hence contributes a fixed fraction to the total. 

From the scaling form of the energy power spectrum ( I20l l, 
the energy decay law ilEl and the length scale growth law 
( I24> one can easily show that a causal power spectrum with 
n = 2 behaves as 



EM{k,t) - kH^ 



(25) 



at low k. This growth at large scales should be taken into ac- 
count when calculating the strength of primordially-produced 
fields, and will weaken the strong bounds derived by Durrer 
et al. (1998, 2003), Caprini & Durrer (2001), Caprini et al. 
(2003). 

5. Comparison with 2D fluid turbulence 



The linear relation between / and Re^ implies that there 
should be a linear relation between s and Re^, and hence a 
quadratic relation between Q{t), the energy decay exponent 
defined in Eq. ( II 9> , and Re^. Fig.|3j showing Q and Re^ 
for 40 < i < 80, confirms that this is indeed the case, with 
asymptote at large Rem consistent with Q — 1/2. Further- 
more, it also shows that both Q and Re^ change little over 
the last half of the runs. 

Let us now turn to the decay of magnetic energy which we 
use to show that the characteristic length scale of the field ^ 



There is an important similarity between 3D MHD and 2D 
fluid turbulence: both exhibit inverse cascade behavior asso- 
ciated with the existence of ideal invariants (Biskamp 1993) 
In 2D fluid turbulence the energy shows an inverse cascade 
(Frisch & Sulem 1984), which is similar to the inverse cas- 
cade of magnetic helicity in 3D MHD turbulence. There are 
also some corresponding similarities between the decay laws 
in decaying MHD turbulence and decaying 2D fluid turbu- 
lence. Chasnov (1997) studied 2D velocity fields initialized 
with a Gaussian distribution, whose energy spectrum E{k) is 
a power law at low k and has an exponential cut-off at high 
k, much as we do. 
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From the Navier-Stokes equation for an incompressible 
fluid one can easily show 



5(U / o / 2^ 



dt 



(26) 



dt 

where u; = V x u is the vorticity (which is a scalar in two 
dimensions, i.e. lj = ojz). We define two length scales from 
the enstrophy, {uj'^)/2, and the palinstrophy, {Wu)'^)/2, 



(27) 



in addition to the diffusion length scale Idis = v >^t- Postu- 
lating power-law decays for the energy and enstrophy, 



(28) 



we find that m — n — 1, and that 

n = -2vt{Lj^) / {v?) = -2{his/lf. (29) 

Hence, a power-law decay again requires that I ~ i^/^, and 
in the limit that the length scale of the flow becomes large 
compared with the diffusion scale, the energy is conserved. 
This is analogous to the conservation of magnetic helicity in 
3D MHD. 

However, the Reynolds number, approximately constant 
in our simulations, behaves rather differently in 2D fluid tur- 
bulence. Defining Re = ul/v, with u — (u^)^/^, we find 
Re ^ t(i+2)/2^ at least for large Reynolds numbers. For small 
Reynolds numbers we expect Re to decay to zero through 
dissipation, pointing to the existence of a critical Reynolds 
number Rcc, at which the system can remain for long peri- 
ods, if initialized close enough to the correct value (Chasnov 
1997). 

The energy spectrum of decaying turbulence is also found 
to be self-similar in 2D. Similar arguments to those presented 
in Sectionl^show that self-similarity again implies I ^ t^^^. 
These arguments also show that if the palinstrophy length 
scale Ip is of the same order of magnitude as the enstrophy 
length scale I, and that enstrophy is being dissipated at the 
diffusion scale ^diff , the high-fc exponent of a power-law en- 
ergy spectrum must be between —4 and —5. Intriguingly, the 
results of Chasnov (1997) seem to indicate a power law some- 
what steeper than the k^^ behavior expected in stationary tur- 
bulence Kraichnan & Nagarajan 1967, Batchelor 1969. 

6. Discussion and conclusions 

We have studied the evolution of decaying 3D MHD tur- 
bulence involving maximally helical magnetic fields. For fi- 
nite magnetic diffusivity there emerges an important quantity 
s = (Cdiff/^n)^, where is the helicity scale defined in Eq. 
(I14t . and ^diff is the diffusion scale. We find — vt, where 
V is the RMS velocity, and hence that s oc Re^, the mag- 
netic Reynolds number evaluated using the helicity scale. The 
magnetic field coherence length (which can be equally well 
expressed as the integral, helicity or relative helicity scales) 
goes as ^ ~ t^/^, magnetic heUcity Hm ^ t^^'^ and magnetic 
energy Em ^ t^^^^^^". A corollary is that Reu is constant 
once the system has reached self-similarity. Furthermore, we 
can extrapolate to the limit of very large magnetic Reynolds 



numbers, useful for example in the early Universe, to find H 
constant and Em ~ t^^^^. 

Our model for the scaling laws should be compared with 
that of Biskamp & Muller (1999). The first difference is that 
they assumed that the non-linear term in the evolution equa- 
tion for the magnetic field was the dominant source of energy 
loss for the magnetic field, and that the magnetic field was 
asymptotically the dominant contributor to the total energy 
E, expressed as F = E^ /Em ^ 1- Then we can write 

E^T^E^/^/^, (30) 

where ^ is a length scale of the magnetic field. They then 
found the phenomenological relation F ~ E/H, which, 
when coupled with E ~ H/^ and the conservation of H, 
gives E ^ t^^/^. 

We emphasize that this model is not inconsistent with 
ours. We infer Ey ^ from the relation — vt, and 
hence that F - <-i/2+2s _ {E/H)t^''. The difference be- 
tween Biskamp & Miiller's assumed relation F ~ E/H 
and ours is small at large magnetic Reynolds numbers where 
s 0. Furthermore, both approaches need to assume only 
that the non-linear and dissipative terms in Eq. J21> are not 
sub-dominant (rather than dominant) and it turns out that both 
scale with time in the same way, as EM/t. In our simulations 
dissipation typically accounted for about 60% of the energy 
loss in the period i ~ 40 to t ~ 80, which means that the 
field is not force-free. 

We believe that our model has certain advantages, in that 
the assumptions going into it give more physical insights than 
the phenomenological (and dimensionally incomplete) rela- 
tion T ^ E/H. Our assumptions are that the magnetic power 
spectrum exhibits a self-similar form i20\ . with power-law 
behavior k^^ at high fc, that resistive dissipation occurs pre- 
dominantly at the diffusion scale ^diff, that there is a separate 
helicity scale ^h, from which it follows that 2 < z < 3. 
We also assume that the eddy turn-over time £^/v is t, and 
that the relative helicity is asymptotically constant, i.e. that 
the energy decays as fast as possible as is consistent with he- 
licity conservation. Given these assumptions, it follows that 
the dissipative and non-linear terms in the magnetic energy 
loss equation ( I2U scale the same way, and we can infer that 
^ ^ t^/'^. From this we derive a scaling law for H, finding 
as expected that H is conserved in the limit of large mag- 
netic Reynolds number, and with the assumption of constant 
relative helicity, we obtain that E ^ in the same limit. 

In summary, for both 3D helical MHD and 2D fluids, the 
key to understanding the free decay of turbulence is self- 
similarity, coupled to a separation between the scale of the 
flow and the diffusion scale. 

Note added. Since the first draft of this work appeared as 
astro-ph/02091 19 vl, there have been several developments. 
A new review of magnetic fields in cosmology has appeared 
(Giovannini 2004). Using different scaling arguments, Cam- 
panelli (2004) has rederived our scaling laws for decaying he- 
lical turbulence, and extended the analysis to the non-helical 
case, where the ideal limit has Em oc t^^, Ey cx t^^, and 
(X t^/^. Extensive 3D numerical simulations, including the 
interesting case of Prandtl number larger than unity, have 
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been carried out by Banerjee & Jedamzik (2003, 2004). The 
authors did not seem to be aware of our work or of Campan- 
elli (2004) and did not attempt to analyse the decay of the 
hehcity, kinetic or magnetic energies in its light. However, 
examination of their Figs. 1 and 7 in the later paper show de- 
cay laws close to t^'^-^ and for the helical and non-hehcal 
cases respectively (with Prandtl number unity). They do not 
quote Reynolds numbers for the simulations so an exact com- 
parison cannot be made. There has also been further work 
on effects of primordial fields on the CMB (Kosowsky et al. 
2005, Kahniashvili & Kahniashvili 2005) and on the evolu- 
tion of magnetic fields in the post-recombination era (Sethi & 
Subramanian 2004). 
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